Astronomy & Astrophysics manuscript no. aa0346 


February 2, 2008 


(DOI: will be inserted by hand later) 





Color-Induced Displacement double stars in 

SDSS 



D. Pourbaix 1 ' 2 '*, Z. Ivezic 1 ' 3 , G. R. Knapp 1 , and J. E. Gunn 1 

o . 

1 Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544-1001, 
» , ■ USA; pourbaix,ivezic,gk,jeg@astro.princeton.edu 



> 

On 



Oh: 

6 

H 



2 Institut d'Astronomie et d'Astrophysique, Universite Libre de Bruxelles, CP. 226, 
Boulevard du Triomphe, B-1050 Bruxelles, Belgium; pourbaix@astro.ulb.ac.be 

3 H.N. Russell Fellow, on leave from the University of Washington 



Received date; accepted date 

Abstract. We report the first successful application of the astrometric color-induced dis- 



CO _ 

placement technique (CID, the displacement of the photocenter between different band- 
, passes due to a varyin g contrib ution of differently colored components to the total light), 



Wielen 



1996) for discovering unresolved binary stars. Using the 



originally proposed by 
Sloan Digital Sky Survey (SDSS) Data Release 1 with ~ 2.5 10 6 stars brighter than 21 m in 
the u and g bands, we select 419 candidate binary stars with CID greater than 0.5 arcsec. 



The SDSS colors of the majority of these candidates are consistent with binary systems 
including a white dwarf and any main sequence star with spectral type later than ~K7. The 
astrometric CID method discussed here is complementa ry to the photometric selection of 
binary stars in SDSS discussed by 



Smolcic et al 



2004), but there is considerable overlap 



(15%) between the two samples of selected candidates. This overlap testifies both to the 
physical soundness of both methods, as well as to the astrometric and photometric quality 
of SDSS data. 



1. Introduction 



It is believed that 50% of all stars belong to multiple systems (!HeintJl96 9) Nc\ eilheless. being 
aware that a specific star is a binary is always useful because either one throws it out of the sam- 
ple or updates the model to describe it and begins some follow-up observations! Whether one 
deals with stellar evolution or galactic dynamics, binaries always receive some special consider- 
ations. So, it is important to be able to detect the binary nature of a star at an early stage of an 
investigation by any possible means. 
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Besides spectroscopy, photometry, and interferometry, astrometry coupled to photometry has 
lately emerged as a way of revealing the binary nature of a source I Wielenll996 ). That method re- 
lies upon either a change in the position of the source as its brightness varies (Variability-Induced 
Movers, VIM) or a photometric-band dependence of the position (Color-Induced Displacement, 
CID). Whereas VIM requires several observations along the brightness variation cycle and only 
works when at least one component is variable, it only takes one image in each band to identify 
a CID and can be done for non-variable stars. Despite the number of multi-band photometric 
surveys, none so far has carried sufficiently accurate astrometry in at least two distinct bands. 
This lack of observations has been recently alleviated by the Sloan Digital Sky Su rvey. 

YorketalLoOO . 



The Sloan Digital Sky Survey (SDSS; 



Abazajian et al. 2003 



and references 



therein) is revolutionizing stellar astronomy by providing h omogen eous and deep (r < 22.5) 



12001 : 



Smith et al. 



photometry in five passba nds (u, g, r, i, and z: lFukugi jgjstjilJl 92 6 



2002) accurate to 0.02 mag (Ivezicetal. 2< 




}unn et al. 1998 



Hogg et al. 



Ultimately, up to 10,000 



deg 2 of sky in the Northern Galactic Cap will be surveyed. The survey sky coverage will re- 
sult in photometric measurements for over 100 million stars and a similar number of galaxies. 
Astrometric positi ons are accurate to better than 0. 1 arcsec per coordinate (rms) for point sources 



with r < 20.5™ dPier et al 



2003), and the morphological information from the images allows 
robust star-galaxy separation to r ~ 21. 5 m 1 Luoton et alJ2003 ). 



Using the SDSS data, we report on the first successful identification of Color-Induced 
Displacement (hereafter CID) binaries. The underlying ideas of that method are given in Sect. [2] 
Sect.|5]describes the simulation that allowed us to optimize the screening of the data described i n 
Sect. |3 In Sect. [5] we present our results and compare them with those of 



Smolcic et al 



(2004), 



who have recently used color selection to identify a stellar locus made of white dwarf+M dwarf 
binaries. . 



2. Color Induced Displacement 



For any unresolved double star and any photometric filter, the position of the photocenter lies 
between the two components. If the two components have different colors, the position of the 
photocenter will change with the adopted filter. The color-induced displacement is the change of 
the position of an unresolv ed binary depending on the adopted filter (Fig.Q. 
Though 



Wielen 



( 1996) suggested that variability induced motion and color induced motion 



could reveal the binary nature of an unresolved star, only the first approach has been applied to 



date, in the framework of Hipparcos (Wielen 1996, 



vea star, only 
nll996 ; EsJl997 



during the preparation of the Tycho-2 catalogue ( H02 et al. 2000), Tycho and 2MASS ( Skrutskie 



Pourbaix et alJl2003l) . However 



1997) positions were compared. Source duplicity came up as a satisfactory explanation for most 
of the discrepant positions (S. Urban & V. Makarov, priv. comm.). 

In the cas e of SDSS, the position of any object together with its magnitude are measured 



in five bands dPier et alJ l2003). Even after the chromaticity effects have been accounted for, the 
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Fig. 1. Position of the photocenter in the different SDSS bands. For double stars, the positions are 
aligned with the two stars and their order follows their central wavelength. Measurement error 
prevents the positions from being perfectly superposed/aligned for a single/double star. The true 
position of the star(s) is represented as a five-branch 'star' 



five positions of a single star do not superpose exactly owing to measurement error (left panel 
of Fig. For double stars, u and z photometric bands will yield the two positions which are 
the most separated because of the largest central wavelength difference between these two filters 
(right panel). 

The position of the photocenter follows the peak of the efficiency of the filter. Therefore, the 
photocenters are not only aligned but also ordered by the filter effective wavelengths. The r band 
therefore plays a central role even if the r photocenter does not necessarily lie right at the middle 
of the u and z photocenters. Instead of requiring that the u-, r-, and z -photocenters are aligned 
and well ordered, one can require that the angle measured from the r-position between u and z 
is 180°. 

From now on, we will refer to the (u, z) angle as the angle measured from the r-position be- 
tween u and z . The distance between it and z (noted 1 1 (u, z) || ) will refer to the angular separation 
between the position of the photocenter measured in the u and z band, respectively. 

We do not use the g- and i-positions because the scatter on the g- and i-photocenter is not 
significantly smaller than on the u and z positions, thus leading to a lower signal-to-noise ratio 
(S/N) as far as || [g, i) || is concerned. 
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3. Simulations 

In order to develop expectations for the relevant parameter space, we first carry out some sim- 
ulations with noise properties consistent with th e actual data, to estimate a lower bound of the 



separation we can expect to notice. According to 



Pier et al 



(2003), the standard deviation on the 



relative position, i.e. band-to-band, for objects brighter than r ~ 20 mag is 31 mas in u and 27 
mas in z but the authors did not quote any correlation between the residuals in those two bands. 

A parent population of ~ 2.4 10 6 stars with u < 21 and g < 21 and good photometry (see 
next section for details) was therefore used to update the precision of the residuals and to derive 
the correlation between them. The residuals in right ascension in u and z have a standard devi- 
ation of respectively 36 and 20 mas, with a correlation of 0.16. In declination, these precisions 
are 39 and 22 mas for u and z with a correlation of 0.18. 

For the simulation, the origin of the offset is placed at the position in the r band and is as- 
sumed to lie in the middle of the segment joining the true photocenters in u and z. The distribu- 
tions of the angle and separation between the u- and z -positions are derived from 200 000 model 
positions generated for both photocenters using the above standard deviations and correlations. 
Such distributions for separations of 1", 20 mas and mas are plotted in Fig. [2] The spread in the 
angle increases as the true separation goes to zero. Using a 99.9% confidence level Kolmogorov- 
Smirnov test, we reject the hypothesis that the distributions of (it, z) for single stars and binaries 
with separations below 10 mas are different. It is worth noting that in the latter case, the corre- 
lation between the u- and z -positions prevents the angles from being uniformly distributed over 

0-7T. 

Though the percentage o f binaries is usually thoug ht to be quite high (ranging from 85% 
among OB stars i HeintJl969 ) down to 30% for M stars I Marchal et alJ2003 ) and even 10% for 
F7-K systems with periods less than 10 years I Halbwachs et ali2003L . very few binaries actually 



induce a noticeable shift of the photocenter. While the color induced displacement is larger the 
more different the colors, too large a magnitude difference causes the fainter component to be 
undetectable. Thus only a small fraction of true binaries will show detectable CID. Coupled to 
the spread of the angle due to measurement error, this means that one should not expect a large 
deviation of the distribution with respect to that of a pure single star population, especially if one 
looks at the whole SDSS sample. 



4. Data 

In this study, we use the public SDSS data release 1 (DR1) which contains a few million stars. 
Because errors are larger towards fainter magnitudes, we selected stars to be bright enough to 
have good astrometry (u,g < 21), thus reducing the number of stars to ~ 2.5 10 6 . We further 
impose 




Fig. 2. Model distribution of the angles and separations between the photocenters in u and z 
assuming a true separation of l"(left column), 20 mas (central column) and mas (right column). 
Each line in the contour (top panels) represents a linear 10% increment. 

1. the condition that the displacement between the u- and z -positions is larger than 0.2". This 
criterion is expressed as: 

V(Aa u - Aa z ) 2 + (AS U - AS Z ) 2 > 0.2 

where the four offset quantities are readily available in the SDSS database; 

2. the precision on the magnitudes is better than 0.1 mag in u and r and better than 0.05 mag 
for g, i and z. 

Because astrometry and photometry are strongly tied together, bad photometry is likely to show 
up as poor astrometry anyway so these constraints on the photometric precision are actually 
safeguards for the astrometry as well. These two criteria yield a sample of only 19982 entries, 
~ 0.8% of the previous sample. 

The angles versus the separations between the u- and z-photocenters the distribution of the 
latter are plotted in Fig.|5] With respect to the simulations from Sect.|51 there is a strong excess of 
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Fig. 3. Observed distribution of the angles and separations between the photocenters in u and 
z. The upper-right corner delimited by the dashed lines corresponds to our revised selection 
criterion. 



points with a large displacement but with a angle. This very striking feature is due to asteroids. 
Both the displacement and colors of all these 2 180 objects are consistent with the asteroids 



already identified in SDSS lllvezic et al 



2002; 



Juric et al 



2002). In the case of asteroids, the 



position of the photocenter changes from u to z because of the genuine displacement of the 
object on the sky in between the two exposures. Because of the way the u and z CCDs are 
located on the detector (the scanning order is rizug), the positions in u and z are on the same 
side with respect to the r-positions, thus yielding a null angle. 

In the single star simulation, the distribution of the angle between the u and z -photocenters 
shows a continuous decrease towards it. However the lower panel of Fig.[3]reveals that, once out 
of the asteroid region of the distribution, the number of large angles actually increases, which is 
consistent with the presence of binaries. 

Even though the distribution shows an increase from 0.2 rad up to n and the simulation 
indicates that no single star is likely to cause a displacement larger than 0.3", we conservatively 
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adopt 1.5 rad and 0.5" for the lower bounds of the angle and separation respectively. This final 
cut leads to 419 candidate binaries. There is no noticeable clustering of these objects in chip 
coordinates nor in (a, S) which rules out the possibility of an instrumental (e.g. corrupted CCD 
row or column) or observational (e.g. unnoticed bad seeing) problem which would have caused 
the displacement. 



5. Results 



5.1. WD+MD bridge 



If the components of a binary have different colors, the resulting colors of the system do not all 
match single star ones. Whereas most combinations of stellar colors remain consistent with the 
stellar loc us, any departur e from the latter can be interpreted as the signature of the duplicity of 



the source. 



Smolcic et al 



(2004) used that color selection to identify a second stellar locus in the 
SDSS data, namely a bridge in the u — g vs g — r color-color diagram between white and M 
dwarfs. 

Our 419 candidate binaries are displayed in the same color-color diagram in Fig.|4] Although 
the two techniques are rather orthogonal, our purely astrometric approach reproduces the essence 



of the color-based results after 



Smolcic et al 



(2004). Such a result was expected since those two 
groups of stars have the largest color difference but are of similar absolute magnitudes and are 
therefore the most likely to yield a noticeable displacement between the u- and z-photocenters. 
Only 61 objects match the criteria of brightness and colors imposed by Smolcic et al, i.e. 7% of 
their bridge stars match our astrometric criteria. 



In order to explain the bottom-left part of the diagram, one cannot rely on the 



assumption 



of a un ique color for all the white dwarfs, especially in u — g and g — r. Following Har ris et al 
(2003) we instead adopt three different model WDs to cover the range of u — g, yet keeping 
r — i = i — z = 0. The coordinates of the three WDs in the (u — g,g — r) plane are respe ctively 



Smolcic et al 



(2004) lies 



(-0.5, -0.6), (0.2, -0.35), (0.5, -0.25). The WD colors adopted by 
between our second and third models. The flux ratio between the white dwarf and the M dwarf 
is assumed to range between 0.01 and 100.0. The three resulting tracks are plotted as thick lines 
in Fig. H 

Since at least one of our extended bridges passes over the quasar region of the (u — g,g — r) 
diagram, could it be that some of our candidates binaries are actually QSOs? From the i — 
z of these points (Fig . |5}, one can conclude that it is seldom (if ever) the case (according to 



Richar ds et al. 



(2001), QSOs have i 



0). All the questionable points are hence consistent 



with a white dwarf + M dwarf pair. 
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u-g 

Fig. 4. Color-color diagram of the putative binaries (triangles) superposed over the original par- 
ent population of 284503 stars (contours). The thick/thin lines represent systems with a M 
dwarf/K7V component. The short thick line close to the center corresponds to A0V+K5III sys- 
tems. Triangles with a circle around have weird colors that could be the cause of the displacement. 



5.2. Additional models 

Owing to the narrow convergence of the three tracks at the M dwarf end, only 50% of the points 
are bounded by the previous models. What other combinations of stars would produce system 
colors consistent with the data? As we have already stated, the colors have to be rather different 
and yet the magnitudes to be rather similar. Assuming a WD companion, the bluer the main 
sequence component, the larger the magnitude difference and the more similar the colors. 

If the M dwarf is replaced with, say, a K7 main-sequence star (assum ing 2.32, 1.01, 0.32, and 
0.15 for u — g, g — r, r — i, and i — z respectively ba sed on the list after 



convolved with the SDSS filters iFukugita et alii 9961) X the flux ratio in the r band becomes a bit 



Gunn&Strvker(1983) 



more constrained: 2.5F Ti k < F r yvD < 250F ri K- The different tracks leading to the same three 
model white dwarfs as in Sect. 15. H are plotted as thin lines in Fig.0] Because of the flux ratio 
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constraint, these tracks do not go down to the WD region but that part of the diagram is already 
covered by the WD+MD model anyway. 

Almost 90% of the candidate binaries can thus be explained with a model involving a WD 
and any main-sequence star redder than K7. What about the remaining 10% whose (u — g, g — r) 
cannot that easily be explained with any of the previous scenarios? Even though several of these 
points are located around (1.1, 0.4) in the (u — g, g — r) diagram, i.e. they are consistent with 
the peak of the stellar locus and hence with a possible contamination by weird single stars, that 
region is totally depleted in the (u — g,i — z) diagram (Fig.|5}- 

Besides the binaries with a WD component, what other systems could exhibit large displace- 
ments between the u and z filters? A AO main-sequence and a giant K5 have rather similar 
magnitudes, yet rather different colors (almost as much as the WD+MD pairs). The flux ratio 
in the r band between the two components is constrained to the range 0.63-6.3. For the K5 gi- 
ant (resp. AOV), we adopt the colors 3.40 (1. 01), 1.35 (-0.21), 0.56 (-0.14), and 0.33 (-0.1 1) for 



g — r, r — i, and i — z respectively ( Gunn & Stryker 



1983 



Fukugita et al. 



1996). The 



corresponding track is displayed in both Fig.[4]and[5](left panel) as a thick line. It is worth noting 
that this track is consistent with the stellar locus in both diagrams thus confirming that the CID 
binaries would not all be identified as outliers in a color-color diagram. 

Eve n though 2 MASS does not go as faint as SDSS, the overlap between the two surveys 



(Finlator et al. 2000) is nevertheless large enough to obtain combined color-color diagrams for 



241 of our putative double stars (~ 60%). Among the combinations of filters, the u — g vs 
J — Kg exhibits the largest departure of these binaries from the stellar locus (right panel of 
Fig-E}. Instead of a bridge as in Fig.|4] the binary locus appears as a narrow horizontal band in th 
combined 2MASS SDSS diagram. Unlike the SDSS colors, the 2MASS color allows us to rule 
out the possibility for the M stars to be giant rather than dwarf. Note also that there are still CID 
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Fig. 6. Left panel: Ratio of z)\\ over the actual angular separation of the two components. 
The thick (resp. thin) lines represent systems with a M dwarf (resp. K7V) component. The short 
thick line in the upper right corresponds to A0V+K5III systems. Central panel: Distribution of 
the g — r color of the parent population (thin line) and the binaries (thick line). Right panel: 
Color-magnitude of the putative binaries 



binaries whose colors are very consistent with the main stellar locus even when 2MASS bands 
are used as well. 

5.3. Contamination 

About 15 stars have weird colors (triangles with a circle around in Figs. |4] and [5}, especially in 
i — z, which is below -0.5 for six of them. In such cases, the noticed displacement would be 
caused by the wrong astro me trie transformation induced by the unusual colors rather than by a 
true displacement of the photocenter. However spectra available for some of these troublesome 
points confirm the duplicity. So, in the worst case, the contamination rate does not exceed 3%. 

5.4. Angular separation 

The likelihood of detecting a CID binary clearly depends on the difference in colors of the com- 
ponents but it also depends on the actual angular separation of the two stars. From the only cri- 
terion adopted so far, namely || (u, z)\\ > 0.5", what can we infer about that angular separation? 
The ratio of || (u, z)\\ to that separation is plotted in the left panel of Fig. [6] 

All the white dwarf systems show a maximum displacement above g— r = 0.5, corresponding 
to 60%, so any separation larger than 0.85" fulfills the criterion. That is consistent with the 
distribution plotted in the central panel of Fig [6] It shows that whereas the parent distribution 
peaks below g — r = 0.5, the binary distribution peaks well above that value. 

Whereas the fractional displacement starts decreasing above g — r = 0.7, the binary distri- 
bution keeps growing up to g — r ~ 1.3 where the displacement has already decreased to 70%. 
This means that the decrease of the fractional displacement is compensated by the actual angu- 
lar separation of the components. Assuming a constant linear separation, these systems should 
thus be closer. A color-magnitude diagram (right panel of Fig.|6) gives credit to that explana- 
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tion. Indeed, the redder, the brighter while the absolute magnitude goes up thus meaning that the 
reddest objects are on average closer to us. 



6. Conclusions 



Wielen ( 1996) as a way of detecting bina- 



The color induced displacement method described by 
ries has been successfully applied to the first public release of the SDSS data. We identify about 
400 systems whose changes in position are essentially consistent with a white dwarf coupled to 
a lower end (later than ^K7) main-sequence star. We therefore expect ~ 2 000 CID binaries at 
the completion of the SDSS observation campaign. 

This ident i ficatio n of binaries is an independent confirmation of the color based results of 



Smolcic et al 



( 2004). However, whereas they had a lower bound on g — r of 0.3, the astrometric 
criterion allows us to identify candidate binaries down to g — r = —0.4. On the other hand, 
color selection they utilized is more sensitive to binaries with angular separations smaller than 
the sensitivity of the CID method. 

Though the approach has proven to give results, its efficiency is extremely low. Whereas 



Marchal et al 



(2003) quote at least 30% of binaries among M stars (the percentage grows with 
the mass of the star along the main-sequence), only 0.02% are detected through their CID effect. 
It is noteworthy that with such a low fraction, the CID binaries do not affect the overall SDSS 
astrometric precision. 

Because of their much better astrometric precision (typically a few /xas), space-based as- 
trometry missions like SIM and Gaia will eventually supersede the SDSS results presented here. 
According to a Gaia preparatory study ( lArenou & Jordi200ll) . the latter could, for instance, de- 
tect a M0 companion to a GO dwarf star at a 3er level at a separation as low as 2.3 mas. In terms 
of separations, this is ~ 200 times better than the sensitivity of the CID method applied to the 
SDSS data. 
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